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Abstract 

We present the detail of a newly developed relativistic transport model (ART 
1.0) for high energy heavy-ion collisions. Using this model, we first study the 
general collision dynamics between heavy ions at the AGS energies. We then 
show that in central collisions there exists a large volume of sufficiently long- 
lived superdense hadronic matter whose local baryon and energy densities 
exceed the critical densities for the hadronic matter to quark-gluon plasma 
transition. The size and lifetime of this matter are found to depend strongly 
on the equation of state. We also investigate the degree and time scale of 
thermalization as well as the radial flow during the expansion of the super- 
dense hadronic matter. The flow velocity profile and the temperature of the 
hadronic matter at freeze-out are extracted. The transverse momentum and 
rapidity distributions of protons, pions and kaons calculated with and without 
the mean field are compared with each other and also with the preliminary 
data from the E866/E802 collaboration to search for experimental observables 
that are sensitive to the equation of state. It is found that these inclusive, 
single particle observables depend weakly on the equation of state. The dif- 
ference between results obtained with and without the nuclear mean field is 
only about 20%. The baryon transverse collective fiow in the reaction plane 
is also analyzed. It is shown that both the fiow parameter and the strength of 



the "bounce-ofF" effect are very sensitive to tlie equation of state. In partic- 
ular, a soft equation of state with a compressibihty of 200 MeV results in an 
increase of the flow parameter by a factor of 2.5 compared to the cascade case 
without the mean field. This large effect makes it possible to distinguish the 
predictions from different theoretical models and to detect the signatures of 
the quark-gluon plasma which is expected to significantly soften the equation 
of state. 

PACS number (s): 25.75.+r 
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I. MOTIVATION 



The main purpose of relativistic heavy-ion colhsions at future RHIC and LHC energies 
is to create in the laboratory a quark-gluon plasma (QGP) and to study its properties. On 
the other hand, for heavy-ion collisions at Brookhaven's AGS the emphasis has been in 
studying the properties of hot dense hadronic matter and the collision dynamics. However, 
the possibility of a phase transition to the QGP in central heavy-ion collisions at AGS 
energies has also been mentioned recently. 

At AGS, a vast body of data have already been collected and analysed during past few- 
years Comparisons of these data with the predictions from theoretical models, such as 
the RQMD 0, the ARC [|], and the QGSM [Q, have revealed much interesting physics. In 
particular, a picture of nearly complete stopping of baryons in central heavy-ion collisions 
has emerged from these studies. Furthermore, it has been shown that baryon and energy 
densities up to 10 times that in normal nuclei have been reached in these collisions, leading 
thus to the suggestion that in central collisions at AGS energies the QGP may have already 
been formed |Q. It has been further shown, based on a relativistic nucleation theory, that 
perhaps one in every 10^ or 10'^ events undergoes the phase transition. However, the above 
suggestion and many other exotic phenomena depend crucially on the maximum energy and 
baryon densities reached in the collision. In this paper, we will carry out a detailed study 
on the formation of a superdense hadronic matter in heavy ion collisions at AGS energies. 

With the available Au beam at pbeami A=11.Q GeV/c at AGS, more systematic study of 
heavy- ion collisions are being carried out by several collaboration (e.g. In addition, 

systematic and exclusive measurements of the energy and mass dependence of particle pro- 
duction, correlations and collective flow effects from Bevalac energies up to AGS energies 
will soon be carried out by the EOS collaboration 0. Besides studying possible new physical 
phenomena, these experiments will provide a broader base of data so that more stringent 
tests of theoretical models can be made. Also, comparisons of the experimental data with 
reliable model predictions will allow us to learn about the properties of hot dense hadronic 
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matter formed in these collisions and help identify new physical phenomena. Stimulated by 
the success of the RQMD and the ARC in describing many available experimental data at 
AGS energies, we have recently developed a relativistic transport (ART) model for heavy-ion 
collisions at these energies. A brief report of the model has been given in refs. PJTO[], here 
we will present the detail of the model (ART 1.0) and study several important aspects and 
issues of heavy-ion collisions at AGS energies. 

More specifically, we discuss in section II the detail of the model ART 1.0 and its inputs. 
In section III, we first study the general heavy ion collision dynamics at AGS energies. 
We then discuss the formation of the superdense hadronic matter and its properties as 
well as its dependence on the equation of state. The degree of thermal equilibrium in the 
superdense hadronic matter and the time scale involved are also discussed. In addition, 
main features of the radial flow during the expansion phase, the flow velocity profile and 
temperature of the hadronic matter at freeze-out are also investigated. The model is then 
used to study the transverse momentum and rapidity distributions of protons, pions and 
kaons with and without including the mean field in order to identify observables that are 
sensitive to the equation of state. These results are then compared with each other and 
also with the preliminary data from the E866/E802 collaboration. Finally, we carry out a 
detailed analysis of the baryon transverse flow in the reaction plane and study its dependence 
on the equation of state. The summary and outlook are given in section IV. Those interested 
only in the results may turn directly to section III. 

II. A RELATIVISTIC TRANSPORT MODEL: ART 1.0 

Based on the well-known Boltzmann-Uehling-Uhlenbeck (BUU) model (e.g. ||TTHT3|) for 
intermediate energy heavy-ion collisions, we have recently developed a relativistic transport 
model for heavy ion collisions at AGS energies P,p!0| by including more baryon and meson 
resonances and their interactions. The BUU model has been very successful in studying 
heavy-ion collisions at beam energies lower than about 3 GeV/nucleon. We refer the reader 
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to the review by Bertsch and Das Gupta for more details of the BUU model |TT . 

In the ART model, we have kept the same philosophy and methods as in the BUU 
model, but added some new physics and numerical techniques in order to simulate heavy- 
ion collisions at higher energies. More specifically, we have included in ART 1.0 the following 
baryons: N, A(1232), A^*(1440), A^*(1535), A, E, and mesons: vr, p, u, r], K, as well as 
their explicit isospin degrees of freedom. Although antiparticles and heavier mesons have so 
far not been included, their production can be studied perturbatively in the present version 
of the model. We plan to include explicitly these particles in an upgraded version of the 
model. Both elastic and inelastic collisions among most of these particles are included as 
best as we can by using as inputs the experimental data from hadron-hadron collisions. 
In accordance with this philosophy, almost all parameterized cross sections and angular 
distributions that have been used in the BUU model are replaced by empirical expressions 
based on the double-logarithmic interpolations of the experimental data . However, more 



than 200 reaction channels are listed in the CERN data book [|1^] for nucleon-nucleon and 
pion-nucleon collisions. We certainly have not fully incorporated all these channels. Instead, 
most inelastic hadron-hadron collisions are modeled through the formation of resonances. 
The advantage of this approximation is that the finite lifetime of these resonances takes into 
account partially the effects of the finite formation time for produced secondaries. In the 
following, we present details of treating various hadron-hadron interactions. 

A. Inelastic baryon-baryon interactions 

First of all, we have included in the model the following reactions, 

iVAT ^ A^A, A^A^*(1440), A^A^*(1535), (1) 

NN ^AA, AA^*(1440), (2) 

NN NNp, NNu, AAvr, (3) 

A^A^ ^ AAp, (4) 
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NA^ NN* (lUO), NN* (1535), (5) 
AA ^ A^A^*(1440), ArAr*(1535), (6) 
AA^*(1440) ^ A^A^*(1535), (7) 

and those producing kaons as shown later in Eqs. (|68| ) - ([70|) . The decomposition of the total 
inelastic nucleon-nucleon cross section into different channels shown in the above certainly 
involves some uncertainties and assumptions. Moreover, most of the available data are only 
for the pp reaction. For the np reaction we use as much data as possible. Otherwise, we 
just assume that the cross sections are the same as in the pp reaction. In the following, we 
discuss the methods we have introduced for determining these cross sections. 

As in refs. [12,16,1^, we use the parameterization introduced by VerWest et al. |15] for 



the cross sections of producing a single A or A^*(1440) resonance in processes shown in Eq. 
(1). In terms of the channel isospin cross section cTjj with i{f) = or 1 being the initial and 
final isospins, the A and A^*(1440) production cross sections are given by 



n + A++ 


, 1 
j = o-io + -an. 


(8) 


P + A+) 
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(10) 
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1 1 
= 2^11 + 4^10' 
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p + N*+ 
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p + N*^ 
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) = ^Oroi, 


(13) 


n + N*^ 


) = ^^01. 


(14) 



In the above, aif{y/s) is defined as 

^^^^^^ = ^"^^^ {s*-miy + mlT^^ ^'^^ 
where the parameters a, (3, mo and F are listed in table |. Other quantities in the param- 
eterization are defined as in ref. [|T^. The cross section for the production of a A^*(1535) 
resonance is estimated from the empirical rj production cross section according to [|18| 



a{NN ^ NN* (1535)) ^2a{NN ^ NNri), (16) 

0.102(v^- 2.424) , 
a{pn pnri) ^ 3a{pp pprj) = ^ _ 2 424)^ 

Cross sections for double resonance production in processes shown in Eq. (2) are es- 
timated by subtracting from the measured inclusive 2tt production cross section in the 
nucleon-nucleon collision |14] the contribution from A^A^ — > NNp and the 2n decay of the 



A^*(1440) resonance in the A^A^ A^A^*(1440) reaction. For example, we have for a pp 
collision, 

a{pp^ ^ AA+ ^ AA^*(1440)) 

Ia^a Ia^n* 

= Aa{pp ppn^n') - 2a{pp ppp^) - 0.4cr(pp A^A^*(1440). (18) 

We show in Fig. 1 the experimental pp inelastic cross section (solid line) together with 
its decomposition into contributions from the production of various resonances. It is seen 
that for ^/s > 3.0 GeV the contribution from double resonance production, shown by the 
dotted line, is about one third of the total inelastic cross section. We have assumed that the 
double resonance production cross sections are the same for all channels allowed by charge 
conservation. This approximation is strongly supported by recent theoretical studies based 
on the one-boson-exchange model . 

The cross sections for p and uj production in reactions shown in Eq. (3) are taken 
directly from the experimental data [|l^]. As an example, we show in Fig. 2 the experimental 
cross sections for the reactions pp — ^ ppuj and pp ppp^ . The cross section for quasi-37r 
production A^A^ — > AAvr is taken as the difference between the measured inclusive 37r and 
UJ production cross sections. For a pp collision, we have 

a{pp ^ AAvr) = 3[(T{pp ppn'^n'ii^) + a{j)p — > pm[^T[^T[~)] — a{j>p ppu). (19) 

IaItt 

It is seen from Fig. 1 that quasi-37r production dominates the inelastic cross section in the 
energy region ^/s > 3.5 GeV which corresponds to the cm. energies of colliding nucleon 
pairs in the early stage of heavy-ion collisions at AGS energies. We then attribute the 



difference between the experimental total nucleon-nucleon inelastic cross section and the 
sum of cross sections from reactions shown in Eqs. (1) to (3) as well as the kaon production 
cross sections from reactions shown later in Eq. ( |68D to the quasi-47r production process 
A^A^ — > AAp (thin solid line in Fig. 1). In this way, the total inelasticity of the nucleon- 
nucleon collision is properly treated as only a limited, though large, number of reactions have 
been incorporated. The errors introduced by this approximation for the quasi-47r production 
cross section is small. From Fig. 1, it is seen that the maximum cross section for the process 
A^A^ — i> AAp at {\^)nn = 4.9 GeV, corresponding to the cm. energy of colliding nucleon 
pairs at a beam momentum of 11.6 GeV/c, is about 3.0 mb, which, however, has to be 
compared with the 13, 9 and 5 mb for the final states AAtt, A A + AA^* and A^A + A^A^*, 
respectively. Cross sections for reactions shown in Eqs. (5) to (7) are taken to be the same 
as in the nucleon-nucleon collision having the same center-of-mass energy and total charge. 

Using the rejection method, the masses of A, A^*(1440), A^*(1535) and p are distributed 
according to the modified single or joint Breit-Wigner distributions with momentum- 
dependent widths. For single resonance production, the cross section can be expressed 
in terms of the transition matrix element /iANN'~*N"r 

m'^Qmlf fy^-'^N drrir , n /" '^^ i i2 



a(NN' N"r) = '-^^^^ ■ / ■ / ^ E \MNN'^N"r\', (20) 



SiS 



where Pi{mr) is the modified Breit-Wigner function 

P^('"^) = , . "'^^^T^l y (21) 

In the above, pi and pf are the nucleon momenta in the center-of-mass frame of A^A^' 
and N"r, respectively. The centroid and width of the resonance are denoted by m^o and 
r(mr), respectively. As discussed in detail in refs. [0,^, to obtain an accurate resonance 
mass distribution requires knowledge about the mass dependence of the matrix element. 
For simplicity, we assume that the matrix element is independent of the mass, and the 
mass distribution of the resonance is then given by Pi{mr). The momentum factor pf in 
Pi{mr) plays an important role in collisions at energies near the threshold for resonance 
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production as it suppresses the production of resonances with masses near the maximum 

mass TTljnax = — Tfln- 

For double resonance production A^A^' —>■ rir2, one has 

4m^m^ir,mLQ f r dmridmr2 „ / \ f dQ 



a{NN nr.) = ^ ^^^^^^ " J J -^^^m^u -.2) ■ J \M 



|2 

NN'^rir2 \ j 



(22) 



where P2{'mj.i^ 171^2) is the joint Breit-Wigner function 



™ ^ = Pf ■ "^ri ■ Ti mr2 ■ T2 



with the conditions 



m„ + < mri< \fs - 7/1^2, 

m„ + < mr2< y/s - nin - 171^,. (24) 

In Eq. (p3D , P/ is the momentum of the resonances in their center-of-mass frame. Unhke other 
models for heavy-ion collisions at AGS energies, resonances in our model can be produced 
with a large range of masses. 

For reactions producing a p resonance, one has the similar expression for the production 
cross section. The p mass distribution is 

^^"^'^ ^ im,-o.77y+i.25r^{mpy ^^^^ 

Moreover, in calculating the masses, decay probabilities and formation cross sections of 
baryon resonances, the following momentum-dependent widths are used, 

r(A(1232)) = 0.47gV {^1 + 0.6g^) (GeV), (26) 
r(Ar*(1440)) = 0.20(g/go)' (GeV), (27) 
r(Ar*(1535)) = 0.15(g/go) (GeV), (28) 

where go is the nucleon momentum in the rest frame of the resonance with its mass at 
the centroid. However, constant widths of 0.151 and 0.0084 GeV are used for p and u, 
respectively. 



The reverse reactions should be treated on the same footing. However, in the present 
version of the model, we have limited ourselves to only two-body collisions as indicated 
by the arrows in Eqs. (1-7). In view of the high baryon density reached in heavy-ion col- 
lisions at AGS energies, many-body interactions should become increasingly important as 
the compression gradually increases during the collision. Both the physics and techniques 
of treating many-body interactions are currently of great intreat but unfortunately without 
clear solution. We therefore postpone the inclusion of three-body and possibly many-body 
collisions to a later version of the model. Cross sections for the two-body reverse reactions 
are calculated by the detailed balance relation taking into account the finite widths of the 
resonances. For the absorption of a single resonance N"r NN' we have the cross section 

a{N"r NN') = "^^"^^^^ f f^Y. |A<^._^^'P, (29) 
2gr7Tspf{l + 6nn') J 47r 

where = 4(2) is the spin degeneracy of the resonance r = A(A^*). The factor (1 + 
Snn') takes into account the case of having two identical nucleons in the final state. Using 
N"r-,NN'\'^ = \-M.NN'-*N"r\'^, Eqs. ( pO]) and (^9]) lead to the following detailed balance 
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^(AT'V NN') = ■ \ . ■ ^ 

2grmf.Q (1 + dNN') Pf 

(ry/s—ni]^ rim' \ ^ 

/ ^Pi(m:.) . (30) 

Similarly, for the absorption of two resonances rir2 —>■ NN' we have 

a{rir2 NN ) = — — — / — }^ \Mr,r2^NN'\ ■ (31) 



grl9r2TTs{l + 6NN')Pf J 4:TT f- 

The detailed balance relation in this case then reads as 

mrimr2 1 + 5rir2 Pi 



a(rir2 ^ A^A^') 



Agr^gr^rnl^Qml^Q 1 + 5nn' Pf 
dm'^i dm'^ 



a{NN' r,r2) ■ / ^^P^im'^i,^^)) ■ (32) 



Limits of the mass integration are given by Eq. (0). In the case, Pi{mr) = T{pf5{mr — 
mro)/{2mro) and y/s oo, the above relations reduce to the standard detailed balance 
relations for the production and absorption of resonances with fixed masses 
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a{N"r ^ NN') = \ 4 " ^(NN' ^ N"r), (33) 



grigr2 1 + ^ATAT/ 



a{nr2 ^ A^A^') = — — ■ 7—7^ ■ ^ ■ aiNN' ^ rirs). (34) 



At AGS energies, resonances with large widths can be produced, especially in the early stage 
of the collisions, it is then necessary to use the detailed balance relations of Eqs. ( pOD and 
(^) to treat their annihilation. 

B. Baryon-meson interactions 

One can also separate meson-baryon interactions into the elastic and inelastic parts. For 
the elastic interaction, we treat them via the formation of baryon resonances, 

ttA^^ A,Ar*(1440), Ar*(1535), (35) 
r]N ^ A^*(1535), (36) 

as well as direct reactions, 

7r + A^(A,A^*) ^ 7r + A^(A,A^*), (37) 
p + N{A,N*) ^ p + N{A,N*), (38) 
K+ + N{A,N*) ^ + N{A,N*). (39) 



The standard Breit-Wigner form |^T| of resonance formation in meson-nucleon interac- 
tions can be rewritten as 

■'f> N..,.-,S':r(^.) ' («) 

where go is the meson momentum at the centroid m^o of the resonance mass distribution. The 
mass rrir of the produced baryon resonance is uniquely determined by reaction kinematics. 
The maximum cross sections are given by 

amax(vr+ + p ^ A++) = (T„,ax(7r- + u ^ A' ) =190 mb, (41) 
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fx,„ax(vr° + p ^ A+) = a^ax(vr° + n ^ A") =50 mb, (42) 

amax(7r" + p ^ A°) = a^ax(7r+ + u A+) =30 mb, (43) 

crmax(vr" N*0(1440)) = (T^ax(vr° + n ^ N*0(1440))= 6 mb, (44) 

(TmU^^ + N*+(1440)) = a^ax(7r° + p N*+(1440))= 12 mb, (45) 

crmU^~+P^ N*°(1535)) = a^^^Ti^ + N*0(1535))= 8 mb, (46) 

cTmax(vr+ + r2^N*+(1535)) = a,,ax(7r°+p^N*+(1535))= 16 mb, (47) 

^max(r/ + P^N*+(1535))= a,,ax(r? + r2^N*0(1535))= 74 mb. (48) 

In evaluating CTmax, we have taken into account properly the decay branching ratios of the 
A^*(1440) and A^*(1535) resonances. 

The formation of the three baryon resonances in the reactions shown in Eq. ( |35|) accounts 
almost entirely the vr + iV elastic cross sections at low energies. This is demonstrated in Fig. 
3 and Fig. 4 where the experimental cross sections for the elastic scattering of tt^ + p and 
7r+ + p are compared with the sum of the three baryon resonance formation cross sections. 
It is also seen from the figures that the formation of the A resonance alone is not enough to 
describe even the low energy part of tt + scattering. At higher energies, the elastic cross 
section is mainly due to the formation of higher resonances which are not included in the 
present version of the model. We therefore attribute the difference between the experimental 
elastic cross section and the contribution from the three lowest baryon resonances to the 
direct process it + N —* tt + N. The contribution from direct it + N scattering is shown by 
the dotted lines in Fig. 3 and Fig. 4. 

The experimentally unknown cross sections, such as those for tt^ + A^, 7r + A(A^*) and p + 
A^(A, A^*) reactions, are calculated by using a resonance model that includes heavier baryon 
resonances with masses up to about 2.0 GeV. Neglecting interferences between resonances, 
one has 

^ ^ ■ P Y(25Af + l)(25B + l)(v^-m«)2 + 0.25r2j(totoO' ^ ^ 

The pre-factor 1.3 is introduced so that the high energy part of tt"*" +p data can be fitted and 
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is mainly due to the neglect of interferences. The summation is over all baryon resonances 
with masses up to 2 GeV. The total width r^(totaZ) and partial widths Tji{M + B) of 
heavier A^* and A resonances used in the summation are listed in Table |I| and Table fT\\ , 
respectively. For illustrations, we show in Fig. 5 and Fig. 6 the calculated 7r° + p{n) and 
7r+ + A^(A°) elastic cross sections. For the + A^(A, A^*) scattering a constant cross 
section of 10 mb is used according to the available data. 

The decay of resonances, such as A(A^*(1440), A^*(1535)) ^ vr + A^, N*{1U0) 2n + N, 
A^*(1535) ^ T] + N, p 27T, and u 3it, during each time step dt are treated by the Monte 
Carlo method. The decay probability of a resonance is calculated via 



where 7 = E/m^ is the Lorentz factor associated with the moving resonance. For the 
decay of A^*(1440) and A^*(1535) resonances the actual final state is chosen according to 
the corresponding branching ratios. Namely, i?(A^*(1440) vr + A^) = 0.6, _B(A^*(1440) 
27r + A^) = 0.4 for m(A^*(1440)) > 2m„ + m^, otherwise 5 (A^* (1440) vr + A^) = 1. For the 
A^*(1535) resonance we use 5(A^*(1535) ^ vr + A^) ^ 5(A^*(1535) ^ri + N) = 0.5. Finally, 
A and S are allowed to decay into a nucleon and a pion only at the very end of the reaction. 
It is worthy to mention that the Bose- Einstein enhancement factor (1 + for a pion in 
the final state has not been included in meson+nucleon collisions and decays of resonances. 
This is certainly another aspect to be improved in a later version of the model in view of 
the high pion densities reached in heavy-ion collisions at AGS energies. 

The n + N inelastic collision starts at {^/s)^^N ~ 1.2 GeV due to the ttttN final state. As 
the energy increases the number of final state quickly increases to the order of 200. Most 
inelastic reactions lead to pion and kaon production. The final states can also consist of 
baryon and meson resonances. To treat the inelastic reactions in a numerically tractable 
way, we model them through the production of A, p and uj resonances, i.e.. 



Pdecay = 1 " exp[-dt ■ T^TTlr)/ (7 ■ k)], 



(50) 



TT + A^ ^ A + vr. 



(51) 



TT + A^ ^ A + p. 



(52) 
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Ti + N A + u, (53) 

and those producing kaons shown later in Eq. (|7^). We have thus included only final states 
with at most 4 quasipions. Again, the above decomposition of the cross sections involves 
some uncertainties and approximations. However, as we will discuss later our final results 
are not sensitive to them due to the large number of final-state interactions. 

For two and three pion production in vr + collisions, we assume that they are mainly 
through the production of Avr and Ap, respectively. For example, for the Avr final state in 
the TT"*" + p interaction, we use 



(j(7r"'"j9 — s> ^ Avr) = 0"(7r"'"j9 pn^ 



+ 2(T(7r+p^ A++7r°), (54) 



and for the Ap final state, we use 



aiji^p — > ^ Ap) = 3cr(7r^p pn^Ti^Ti ) 

lAlp 

+ 3a{7T^p pn^p^) 
+ 2a{7T+p A++p°) 

+ a{7i^p ^ puj). (55) 

The difference between the experimental total vr + inelastic cross section and the cross 
sections for producing An, Ap and kaons are then attributed to the production of Au, i.e., 

a(7r+p ^ A+c.) = a:Sj;(7r+p) - a(7r+p ^ 

IaIw 

- a{7T+p ^ ^ Ap) - cx{7r+p K+X). (56) 

lAlp 

As an example, the decomposition of the inelastic 7r+ + p cross section is shown in Fig. 
7. It is seen that the cross sections for having 2, 3 and 4 quasipions are comparable at 
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•y/s f« 2.5 GeV, but at higher energies the Att channel dominates. As we shall see later, the 
distribution of meson+baryon center of mass energies peaks at about ^/s Ri 1.5 GeV in the 
reaction of Au+Au at Aeam/^= 11-6 GeV/c and falls off very quickly towards ^/s a; 2.5 
GeV. It thus makes us feel rather confident that the errors introduced in the quasi-47r channel 
in our previous decomposition is very small. The two-body reverse reactions of the tt + N 
inelastic collisions are similarly treated as for the two-body reverse reactions of the inelastic 
baryon-baryon collisions. 



C. Meson- meson interactions 

We model pion-pion elastic collisions through the formation of a p meson, i.e., tt + tt p, 
and the direct process n + n — > n + n. The latter takes into account the case when the 
quantum numbers of colliding pions forbid the formation of a Pi i meson p. 

In terms of the channel isospin cross sections with /=0, 1, and 2, one finds using 
the Glebsch-Gordan coefficients the following isospin decomposition of tttt scattering cross 
sections 

(57) 
(58) 
(59) 
(60) 
(61) 

Because of the symmetrization of two-pion states, only even partial waves are allowed for 
(7° and (7^ and odd partial waves for a^. Low energy tttt scattering are therefore described 
by the three partial- wave cross sections (Tq, ctq and crj. The formation of a p resonance has 
thus the probability 

Ppi^^) = l-^y (62) 
z a[7nr ) 
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= a{7i n ) — a'^, 


(7(7r+7r°) 


= (7(7r-7r°) = ^a^ + ^a' 




= 1^2 ^ 1 1 ^ 1 

6 2 3 ' 
1 2 


(7(7rV) 






- vrV) = + la'. 



while the direct vrvr vrvr reaction has the probabihty of 1 — Pp(7r7r). 

The partial wave cross sections o"£ are related to the phase shift (5£ via 

ai = ^{2L + l)szn%, (63) 

where q is the momentum of the pion in the center-of-mass of the two pions. For 6^ we use 
the parameterization of ref. 



'°-'^"A^^^)- '''' 

5l= -0.12q/m^, (66) 

where Tp = 0.095q[q/mT,/{l + (g/0.77)^)]^. The isospin-averaged cross section for vrvr scat- 
tering 

ainn) = + ^al + (67) 

is shown in the upper window of Fig. 8 as a function of the center-of-mass energy ^/s. Since 
the densities of p and u are rather small at AGS energies, we have neglected in the present 
version of the model the elastic scattering of vr + p, tt + and p + uj. 

The inelastic collisions among mesons are modeled through the reaction, MM — > KK. 
The cross section for this reaction is highly uncertain, here we use the cross section calculated 
from the i^*-exchange model of ref. E^. The lower window of Fig. 8 shows the isospin- 



averaged cross section for the reaction tttt KK. Calculated cross sections in ref. ||23| for 
the reactions pp KK and up KK have values about 0.3 mb, except at energies very 
close to the threshold, which we shall use in our model. Furthermore, we assume that for 
the reactions ttuj KK and pu KK, the cross sections have similar values. 

D. production in baryon-baryon and baryon-meson interactions 



We are at present mainly interested in K~^ production. The following kaon production 
channels in baryon-baryon collisions are included 
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NN ^ NA{E)K, AA{E)K, (68) 
NR^ NA{i:)K, AA(S)K, (69) 
RR NA{J:)K, AA(S)Js:, (70) 

where R denotes a A, A^*(1440) or A^*(1535). We use the approximation that the kaon 
production cross sections in reactions induced by resonances are the same as in nucleon- 
nucleon colhsions at the same center-of-mass energy which are mainly taken from the data 



compilation of refs. |14.24]. Since the experimental data on kaon production are very limited, 



we use only the isospin-averaged cross sections. Based on the one-pion exchange model these 



cross sections can be expressed in terms of the experimentally known ones [p5 | 



a{NN ^ NAK+) ^ a{NN AAK+) ^ ^a{pp pAK+), (71) 
a{NN ^ a{NN ^ A^K+) ^ ^[a{pp pTPK+) + a{pp pS+ir°)]. (72) 

The threshold energies are 2.56, 2.74, 2.63 and 2.77 GeV for the four final states 
NAK, AAK, NUK and ASK, respectively. The actual final state of a given collision is 
determined using the Monte Carlo method according to the relative ratios among these 
cross sections. 

The kaon momentum distribution from a baryon-baryon interaction was parameterized 
by Randrup and Ko according to a modified phase space |2^, i.e.. 



E dMV~s) E 12 f p \ ( p V ^ 

Pmax \ Pmax / \Pmax / 

where Pmax is the maximum kaon momentum given by 

Pmax = ^\J[s - {ms + my + itlkYWs - {tub +mY - mKY]/s. (74) 

The angular distribution is taken as isotropic in the center-of-mass frame of colliding baryons. 
In meson-baryon interactions, kaons are produced through 

7r + Ar(A,Ar*) ^ A(S) (75) 
p + N{A,N*)^ A{T) + K, (76) 
uj + N{A,N*) ^ A{j:) + K. (77) 

17 



For pion-nucleon collisions, the isospin-averaged cross sections for kaon production can 



be expressed in terms of the available data p7| 



ainN ^ AK+) ^^(r{7i+n^ AK+), (78) 
a{7cN ^K+) ^ ^[cT{7i-p ^-K+) + a{7i-p T.^K^) + (x{7f+p ^+K+)]. (79) 

Since there are no data available for kaon production in resonance induced reactions (e.g., 
p{uj) + A^(A, A^*) and vr + A(A^*)), they are taken for simplicity to be the same as in the 
vr + collision at the same center-of-mass energy. This is probably a resonable, minimal 
assumption one can make. However, it is interesting to note that theoretical efforts have 
recently been made to calculate these unknown cross sections (e.g., [^), and we plan to 
incorporate these results in future versions of the ART model. 

The experimental information about the angular distribution of the outgoing kaon in a 



meson-baryon collision is rather sparse but shows a complicated structure [0. Here, we 
use an isotropic distribution in the meson-baryon center of mass frame. A test using a 
forward-backward peaked angular distribution in the baryon- meson cm. frame shows little 
change in the final kaon distribution in heavy-ion collisions as a result of kaon final-state 
interactions. 

Finally, we would like to point out that in reactions RN NY K and RR — > NYK, 
which are dominated by one-pion exchange, the pion can be on-shell. The latter contribution 
is equivalent to a two-step process, i.e., a resonance decaying into a physical pion and a 
nucleon, and the subsequent production of a kaon from the pion-baryon interaction. To avoid 



double counting, one can follow the approach of Ref. p8[ by including only the reactions 
RN — > NYK and RR NYK but not the reaction tcB —* NYK. This requires, however, a 
model for the reactions -RA^ — ^ NYK and RR —>■ NYK. In the present paper, we use instead 
the assumption that the cross sections for the reactions -RA^ —>■ NYK and RR —>■ NYK 
due to off-shell pions are the same as in A^A^ NYK reaction at the same center-of-mass 
energy, which invloves only an off-shell pion. Then, we should include both the reactions 
RN NYK and RR NYK and the reactions tcB NYK in the transport model. 



E. Hadron momentum distributions 

In this section, we discuss the method we use to determine the final particle momenta 
in a hadron-hadron collision, which does not involve kaons. In principle, to determine the 
momentum of a particle in a multiparticle production process, one needs to know the interac- 
tion matrix element which is mo del- dependent. For simplicity, we determine the momentum 
distribution of the produced particles according to the phase space. Besides imposing the 
energy-momentum conservation, we also make use of the systematics observed in multipar- 
ticle production processes. These systematics are often given by empirical formula that are 
fitted to the inclusive data. 

An important feature of energetic hadron-hadron collisions is the "leading" particle be- 
haviour. In our model this is ensured by requiring the outgoing baryons to have the same 
or similar identities as the incident ones so that their longitudinal directions are retained 
when performing the momentum transformation from the baryon-baryon cm. frame to the 
nucleus-nucleus cm. frame. Moreover, baryons in the final state of an energetic hadron- 
hadron collision have a typical forward-backward peaked angular distribution and a soft 
transverse momentum distribution. For baryon-baryon collisions with ^/s < 3.0 GeV where 
single-resonance production dominates, we use the angular distribution obtained from fitting 
the pp NA data 0], i.e., 

-^ = biis) + 3h{s)cos'e, (80) 

where = 0.5 and bs{s) = for y/s < 2.14 GeV, and 

= 29.03 - 23. 75v^ + 4.87s, (81) 
63(5) = -30.33 + 25.53v^- 5.30s, (82) 

for 2.14 < < 2.4 GeV, while for 2.4 < < 3.0 GeV, bi{s) = 0.06,63(5) = 0.4. This 
distribution is shown in Fig. 9, and it is seen that the forward-backward peaked angular 
distribution becomes almost energy independent for ^/s > 2.3 GeV. This tendency is consis- 
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tent with the systematics found in pp colhsions at high energies where the inclusive proton 



with X* being the scaled longitudinal momentum in the cm. frame, i.e., x* = 2j9*/y^. The 
distribution has the properties of the naive scaling, pt and p* factorization and the soft, 
energy-independent transverse momentum distribution. 

Here we adopt the same functional from of momentum distribution as given by Eq. (|83D 
for all baryons in collisions with y/s > 3.0 GeV, but slightly adjust the parameters in the 
transverse momentum distribution so that we can reproduce the inclusive proton momentum 
distribution from pp collisions at Pbeam = 15 GeV/c. Although both longitudinal and trans- 
verse momentum distributions are needed for three-body final states, only the transverse 
momentum distribution is required for two-body final states as the longitudinal momen- 
tum can be determined from the momentum conservation. The momentum distribution 
determined by the present model turns out to be very similar to Eq. (|83D , i.e.. 



In the upper window of Fig. 10, the proton transverse momentum distribution from pp col- 
lisions at Pbeam = 15 Gcv/c, which is used to determine the input momentum distribution 
in Eq. (^4|), is shown by the solid squares. The solid line is a plot of the proton transverse 
momentum distribution from Eq. (^). The quality of the fit is reasonably good although 
there is still some room at low transverse momenta for further improvement. The longi- 
tudinal momentum distribution has a very similar behaviour. With the above transverse 
momentum distribution for baryons and employing the assumption that all resonances decay 
isotropically in their rest frames we obtain the transverse momentum distribution of vr" as 
shown in the lower window of Fig. 10. The tt~ transverse momentum distribution is then 
seen to fit reasonably well the empirical formula. 



spectra can be well described by |^ 




(83) 




(84) 



d^p* 




(85) 
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of Ranft el al p9| for high energy pp colhsions. For comparisons, we have also shown in 
Fig. 10 the generated proton and tt^ transverse momentum distributions from pp collisions 
at Pbeam = 3.0 GcV/c. They, of course, do not agree with the scaling formula valid at high 
energies. 

The above procedure makes us feel confident about our model for high energy heavy-ion 
collisions. We remark here, however, that there is no direct, transparent extrapolation from 
p+Pi p+A to AA collisions. The dynamics and the final momentum distributions of hadrons 
in heavy-ion collisions are rather sensitive to the modeling or prescription of individual 
hadron-hadron scattering, besides the elementary momentum distributions discussed above. 

F. Nuclear equation of state and causality 

The nuclear equation of state describes the response of the nuclear matter to changes 
in excitation energies and densities. Current knowledge on the nuclear equation of state is 
restricted to a narrow region around the ground state. Extending our knowledge about the 
nuclear equation of state to different densities and/or excitation energies has been a major 
goal of nuclear physics. 

Mean-field effects have already been found to be important in heavy-ion collisions at 
medium energies. At higher energies such as at AGS, they have been ignored in ARC, for 
example. We believe that the mean-field potential is not negligible in heavy-ion collisions 
at AGS energies. Although the forward scattering amplitudes of hadron-hadron collisions 
in the high energy limit have been found approximately purely imaginary ||30|, the AGS 



energies may not be high enough for the real part of the scattering amplitude to completely 
vanish. Of course, the form and strength of the mean field in hot dense hadronic matter 
is uncertain and has been a subject of many discussions. Secondly, although in the early 
stage of the collision the kinetic energy is much higher than the potential energy, particles 
are gradually slowed down and the mean field is expected to play an increasingly important 
role as the collision proceeds. In particular, the repulsive mean field in the high density 
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region tends to keep particles from coming close and therefore reduces the maximum energy 
and baryon densities reached in the collision should there be no mean field. Moreover, in 
the expansion phase of the collision mean-field effects become even stronger. It is therefore 
essential to study how the collision dynamics is affected by the nuclear mean field. Of course, 
the most challenging task is to identify the experimentally observable consequences of the 
nuclear mean field, so the information about the nuclear equation of state can be obtained. 

Without much reliable knowledge about the equation of state in hot dense medium, we 
use here the simple, Skyrme-type parameterization for the mean field which has been widely 
used in heavy-ion collisions at and below Bevalac energies, 

U{p) = a^ + h{^Y + Vcouiomh. (86) 
Po Po 

The corresponding energy per nucleon in nuclear matter at zero temperature is given by 

where Ef=37.26 MeV is the Fermi energy. The parameter a is negative, b is positive and a 
is larger than one reflecting the fact that the nucleon-nucleon interaction has a short-range 
repulsive part and a long-range attractive part. By fitting to the ground state properties of 
nuclear matter the three parameters are determined in terms of the compressibility coefficient 
K as 

K + 44 73 

a = -29.81 - 46.90-^--^^^(Me\/), (88) 

5 = 23.45 ^ ^^^^-^^ MeV), (89) 
i^- 166.32^ ^ ' 

K + 44.73 , , 

a = . (90) 

211.05 ^ ^ 

The compressibility coefficient K can also be expressed in terms of a, b, and a, i.e., 

K = 9p(dP/dp)s = 9(^ + a + ba). (91) 

3m 

An important constraint on the nuclear equation of state in hot dense medium is imposed 
by the causality pT| , P^ . For nuclear matter at zero temperature the adiabatic sound velocity 



IS 
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= -i^)s = -llEfi^f' + a^ + bai^n (92) 
m dp m 3 po po po 

The causal requirement < limits the density range applicable for a given parameter 
set. In Fig. 11, we show the adiabatic sound velocity and the energy per nucleon as functions 
of density for the two equations of state corresponding to K=377 MeV (cr =2, stiff) and 
201 MeV (cr = 7/6, soft). It is seen that the stiff and soft equations of state violate the 
causality beyond densities p > 3po and p > 7po, respectively. At finite temperatures the 
critical densities are reduced by less than one density unit |]3^. In the following, we shall 
use the soft equation of state, and the stiff equation of station will only be mentioned for 
the purpose of discussions. 

The mean-field potentials for baryon resonances are likely to be different from that of 
nucleons. Resonance potentials have been studied previously |^ and are still less uncertain 
3^ than the nucleon potential. For example, significant efforts have been devoted to under- 



standing the effects of the A potential on the nuclear matter saturation properties ||3^-|37 



The A optical potential has also been studied in charge-exchange reactions [^, electron 
and 7 induced reactions [^] as well as vr-nucleus scattering |]40| , |4T| . All studies indicate 
that the shape of the A potential is very similar to the nucleon potential, probably a little 
deeper than the nucleon potential. The exact form of the A potential somehow depends on 
the momentum and density. Current theories give a large range of model parameters for 
the A potential, and they often have conflicting features. In view of the larger uncertainties 
associated with resonance potentials, we assume here that they are the same as the nucleon 
potential. 

Since a large number of baryon resonances are produced in heavy-ion collisions at AGS 
energies, the reaction dynamics might be significantly affected by the resonance potentials. 
Relativistic heavy-ion collisions may therefore provide valuable information about the reso- 
nance potentials. However, this can only be possible if other details of the collision dynamics 
are well understood. The present work is thus also useful in this respect. 
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III. APPLICATION OF ART 1.0 TO HEAVY-ION COLLISIONS AT AGS 

ENERGIES 



We first discuss the range of beam energies in wliicli tlie model outlined above is ap- 
plicable. In the energy range of i^beam/^ < 3.0 GeV/A, our model, as many other nuclear 
transport models, is able to describe many aspects of heavy-ion collisions. These include 
the emission of light particles, the creation of pions, kaons and etas as well as the transverse 
collective flow of various hadrons Since the particle production mechanism used in 

the model is exclusively based on the excitation and decay of several low-mass baryon and 
meson resonances, the model is expected to be inadequate once the energy is above about 
Ebeam/A =15 GcV. At higher energies, such as the SPS/CERN energies, heavier resonances 
and more importantly the formation and fragmentation of strings and ropes are expected 
to be the dominant mechanisms for particle production 

The model ART 1.0 is designed intentionally to be most applicable in heavy-ion collisions 
below the AGS energies. As in other hadronic models, our model is naturally limited only to 
interactions among hadrons. If there were a phase transition to the QGP at AGS, hadronic 
models, such as the present one, are still useful as their predictions provide the background 
against which new physical phenomena can be identified. In fact, one of the motivations 
for the present work is to study whether and how conditions for forming the quark-gluon 
plasma can be achieved in heavy-ion collisions at beam energies from 2 GeV/nucleon up to 
AGS energies. Moreover, hadronic models are useful for studying in-medium properties of 
hadrons. Theoretical studies (e. g. |4^-^) have shown that the properties of hadrons, such 
as the mass and decay width, may be modified in hot dense hadronic matter as a result 
of the partial restoration of chiral symmetry. In the following, we shall apply the present 
model to study several general, but important aspects and issues of heavy-ion collisions at 
AGS energies. In the subsequent work we plan to study among many other subjects: (i) the 
effects of changing hadron properties in hot dense medium on heavy-ion collision dynamics 
and, in particular, on the production of exotic particles, (ii) the beam energy dependence 
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of the collision dynamics from 2 GeV/nucleon up to AGS energies which will be carried out 
experimentally by the EOS collaboration, and (iii) the role of quark and gluon degrees of 
freedom at high densities and the study of physics related to the phase transition to the 
QGP. 



A. Reaction dynamics at AGS energies 

To study the reaction dynamics of central heavy-ion collisions at AGS energies, we shall 
first study in this section the reaction rates of various hadron-hadron collisions. To avoid 
the complexity introduced by the mean field, we shall use the cascade mode of ART. In 
particular, we shall consider head-on collisions of Au+Au at pheam/A=^^-Q GeV/c as an 
example for discussions. 

Fig. 12 shows the reaction rates of nucleon-nucleon scattering to the following final 
states NR, RR, AAtt, NNp and NNuj where R denotes baryon resonances A, A^* (1440) or 
A^*(1535). As one would expect from the decomposition of the total inelastic cross sections 
discussed in the previous section, the rate of the quasi-Svr channel is the highest at the 
very beginning of the reaction when t < 1 fm/c, but is closely followed by the quasi-27r 
and quasi-lvr channels. Since particles loses their energies after several collisions the single 
baryon resonance production soon becomes dominant. The collision rates for producing 
p and uj are, however, about an order of magnitude smaller than the quasi-47r channel. 
Most collisions cease after about 10 fm/c except the single-resonance production channel 
A^A^ —>■ NR which lasts a little longer. For comparison, it is worth noting that the maximum 
reaction rate reached in the collision is about an order of magnitude higher than that in 
heavy-ion collisions at Bevalac and SIS/GSI energies |1T2| , ^ . 



The rescattering or absorption of produced baryon resonances and mesons are very im- 
portant in heavy-ion reactions at AGS energies. In Fig. 13 we show the reaction rates for 
nucleon-l- resonance, resonance-|-resonance and meson+meson scattering. It is seen that the 
nucleon+resonance scattering start earlier and dominate during most of the collision time. 
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The meson+meson and resonance+resonance rescattering start later as it takes some time 
to create these secondary particles. After about 10 fm/c the nucleon+resonance and reso- 
nance+resonance scattering are almost over. One also notices that meson+meson scattering 
last longer mainly due to the large tt + tt elastic cross section, although the rate decreases 
towards the later time of the collision. All these features are what one would expect and 
easily understandable. The relatively low rate for baryon resonance+baryon resonance col- 
lisions needs, however, some explanations. The most important reason is that the final 
states in resonance+resonance collisions are limited only to NN, NR' and K~^X in the 
present version of the model, instead of allowing for the production of up to 4 quasipions 
as in nucleon+nucleon collisions. In addition, due to the large spin degeneracy in the initial 
state, the cross sections obtained from the detailed balance relation for the RR — > NN[R') 
reactions are relatively small compared to their reverse reactions. 

In Fig. 14, we show the decay rate of baryon resonances together with their formation rate 
from TT+N collisions. An interesting feature of these rates is that they last a longer time than 
hadron+hadron collisions due to the large meson+baryon cross sections. Another feature 
is that the small difference between the rates of decay and formation of baryon resonances 
makes the apparent lifetimes of these resonances much longer than their lifetimes in free 
space. This is rather clear from the evolution of the pion and baryon resonance multiplicities 
shown in Fig. 15. It is also seen that about half of the baryons are in their excited states at 
about 4 fm/c, which is the instant of the highest compression as we will discuss in the next 
section. The abundance of baryon resonances is dominated by A (1232) resonances as one 
would expect in this energy range. However, it should be stressed that the population of 
A^*(1440) and A^*(1535) resonances is also significant. These higher resonances serve as an 
energy reservoir and thus play an important role in the production of particles, especially 
those having higher energy thresholds [^. It is also interesting to note that the number 
of pion-like particles (free pions + baryon and meson resonances) is about the same as the 
initial number of resonances at the time of maximum compression. 

In Fig. 16, we show the time integrated total number of hadron-hadron collisions as a 
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function of the center-of-mass energy ^/s of the colhding pairs for N+N, N+R, R+R, M+B 
and MM coUisions, where the meson+baryon (M+B) and meson+meson (M+M) colhsions 
include contributions from both baryon and meson resonances. It is seen that N+N colhsions 
start at about ^/s = 1.8 GeV and peak at about ^/s = 2.2 GeV. The peak is mainly due 
to the large number of elastic scattering. The position of the peak shifts to higher energies 
and its hight decreases substantially in N + R and R + R reactions. The reasons are mainly 
the following. First of all, we have assumed that the elastic cross sections for N + R and 
R + R collisions are the same as for N + N collisions at the same center-of-mass energy 
y/s. Since the elastic cross section decreases quickly as the center-of-mass energy increases, 
the heavier masses of the baryon resonances then result in smaller elastic cross sections for 
N + R and R + R collisions. In addition, as we have mentioned earlier, the number of final 
states in N + R and R + R collisions are reduced compared to that in + collisions. 
For meson+baryon reactions, the collision number distribution have two distinct peaks at 
^/s ^1.2 GeV and ~ 2 GeV, respectively. The first peak apparently corresponds to the 
formation of the A(1232) resonance, while the second peak has contributions also from the 
A^* resonance and direct it ~ N collisions. The meson+meson collision number has a peak 
at about the p mass ~ nip =0.77 GeV. 

The above examination on the collision dynamics indicates that the model is working 
well under control, regardless of the question whether some of the assumptions made in the 
model can be improved. This gives us further confidence about the model. 

B. Formation of superdense hadronic matter 

Depending on the nature of the quark-gluon plasma phase transition at high baryon 
densities, the crossover energy between the hadronic and quark-gluon phases may occur 
anywhere between Eiab/A = 2 — 10 GeV pO| , |5T| . It was pointed out recently that the current 
Au+Au reactions at 11.6 GeV/c may have already overshot the transition region [^. It 
is therefore of great interest to critically examine the baryon, meson and energy densities 
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reached in these colhsions. In this section, we study the local baryon, meson and energy 
densities in head-on collisions of Au+Au at Pbeam/^ =11.6 GeV. It will be shown that in 
these collisions a superdense hadronic matter with a density higher than the critical density 
for the phase transition can be formed in a sufficiently large volume and for a sufficiently 
long time. We will also study the dependence of the size and lifetime of the superdense 
hadronic matter on the nuclear mean field. 



The standard test particle method is used to calculate the global density pg in the 
nucleus-nucleus cm. frame on a lattice of size 40 — 40 — 48 with a volume of 1 fm^ for 
each cell. In the test particle method, one replaces the continuous phase space distribution 
function with a finite number of test particles representing individual phase space cells, i.e. 

f{r, P, ^ ^ E ^i^ - r^mip - pM, (93) 

where and p^ are the coordinates and momenta of the test particles. Nt is the number 
of test particles per nucleon, we usually use Nt = 100 for Au+Au collisions. The global 
baryon, meson and energy density function Pg, and Cg are evaluated via 



i=TT,K,rj,ui 



t) = T.J ^,{p' + m'r'h{r,P,t)+Y. J ^s^iP' + ^'y^'f^i-^P^t)- (96) 



In the last equation, i runs over all particles with fixed masses while j runs over all particles 
with variable masses. 

The local baryon, meson and energy densities in each cell are then obtained from pj* = 
Pg/^' ~ hi — ^ghi where 7 is the Lorentz factor of each cell. 



7 = ./(E^0^-(Ep^)VE^^- (97) 

y « i i 

In the above, i runs over all particles in the cell. For most parts of the colliding nuclei, the 
use of local densities can eliminate the trivial compression effects due to the initial Lorentz 
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contraction of the two nuclei. However, in the most central cell where the two nuclei are in 
touch and streaming towards each other, the 7 factor is almost 1. Therefore, we have at the 
very beginning of the collision, ~ = 270P0 where 70 is the center of mass Lorentz factor 
of the colhding nuclei, and for Au+Au at Pteam = 11-6 GeV/c, one has 70 = 2.6. Fortunately, 
this initial high baryon density is only in a very small region at the very beginning and is 
soon destroyed by the violent collisions that follow. In the cascade model calculations, the 
initial high baryon density in the most central cell is simply a matter of presentation. For 
calculations including the density-dependent mean field we have found from comparisons 
with calculations using the scalar density that the initial high baryon density has no effects 
on the collision dynamics and the final observables. 

We show in Fig. 17 the evolution of the local baryon density distribution in the re- 
action plane using the cascade mode of ART 1.0 for the head-on collision of Au-|-Au at 
Pbeam/ A=11.6 GeV/c. The outmost contour is for pi = 0.5po, and the numbers in the figure 
are the densities of the corresponding contours in units of po. The two Lorentz contracted 
nuclei are set in touch at t=0. It is seen that the initial high density region bounded by 
the contour with p = 2po is very small. The two nuclei soon form a compressed disk and 
reach the maximum compression of about 9po at about 4 fm/c. The compressed hadronic 
matter then starts to expand at about 6 fm/c. We see that up to 10 fm/c the longitudinal 
expansion is significantly stronger than the transverse one. This observation seems to in- 
dicate that the isotropically expanding, spherical fireball model widely used for describing 
heavy-ion collisions at AGS energies may need to be improved. It should be mentioned that 
the main features observed here are in general agreement with those extracted from the 
ARC calculations |^. 

The local densities are well defined and are the relevant quantities for discussing the 
interesting physics of the phase transition to a quark-gluon plasma. However, the current 
model is still not fully relativistically covariant. To estimate the minimum compression 
reached during the collision for a given system at a fixed beam energy and impact parameter, 
we show in Fig. 18 the scalar baryon density pi defined as 
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h, ^ /■ d^P p / ^fn /■ d^p dm ^ , , m 



where e is the energy of a baryon. In the hmit that the size of each cell is very small such 
that it can only accommodate one particle the local density becomes equal to the scalar 
density. It is seen that the main features in the evolution of the scalar density are similar to 
those of the local density. It is interesting to note that the initial high density region around 
the most central cell in the local density plot does not appear here as one would expect. 
We have found that the maximum scalar density reached at 4 fm/c is about 6.5po which is 
significantly smaller than the maximum local density reached in the collision. 

The local density of mesons (mostly pions) can be as high as that of nucleons during the 
collision. The evolution of the meson density distribution in the reaction plane is shown in 
Fig. 19. Mesons mainly populate the central region for t < 2 fm/c, and they spread out 
quickly soon after that and finally occupy the whole collision volume. The maximum meson 
density of about 4.5po is reached at about 4 fm/c. As we have discussed earlier, at this time 
about half of final pions are still bound in resonances. It is also interesting to see that mesons 
start to expand earlier and faster than baryons. After about 8 fm/c, the meson sector of 
the system has already been expanding almost isotropically. In principle, the geometry and 
time scale of the expanding hadronic matter studied here can be compared directly with the 
interferometry measurements to further test the model. 

To form a quark-gluon plasma, it is necessary to achieve a high local energy density in a 
sufficiently large volume and for a sufficiently long time such that the initial plasma droplets 
can grow. The current estimate for the critical baryon and energy densities at which the 
QGP may form are about 5po and 2.5 GeV/fm^, respectively From the above studies on 
the baryon and meson densities, we have seen that the superdense hadronic matter occupies 
a large volume and for a rather long time. As to the possible phase transition to the quark- 
gluon plasma, more complete and quantitative information about the local energy density of 
the superdense hadronic matter are needed. Fig. 20 shows the evolution of the local energy 
density in the reaction plane for the head-on collision of Au+Au at pbeam = 11-6 GeV/c. 
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The numbers in the figure are the energy densities in units of GeV/fm^. It is seen that 
a small, thin disk of energy density higher than 2.5 GeV/fm^ has already been formed at 
about 2.0 fm/c. It grows until about 4 fm/c, then starts to decay, and at about 7 fm/c it 
disappears due to expansions. The superdense hadronic matter has thus a lifetime of about 
5 fm/c and a maximum volume of about 200 fm^ at about 4 fm/c. 

The results obtained from the cascade calculations seem to indicate that in head-on col- 
lisions of Au-|-Au at p6eam/^=ll-6 GeV/c, the conditions for forming a quark-gluon plasma 
has indeed been reached. However, this also marks the breakdown of the hadronic models 
and the need to include new degrees of freedom. It is thus a great challenge to actually 
convert this superdense hadronic matter into the quark-gluon plasma. Without introducing 
explicitly the quark and gluon degrees of freedom in the model make this impossible. Nev- 
ertheless, before one can study the properties of the quark-gluon plasma, it is necessary to 
investigate how and to what degree the properties of the superdense hadronic matter may 
be affected by the nuclear mean field, which is an important piece of physics neglected in 
cascade calculations. In the following, we shall compare results obtained with and without 
the mean field. 

First of all, we have found that the general features in the evolution of the baryon, meson 
and energy density distributions are very similar in the two calculations. However, important 
quantitative differences exist. Fig. 21 shows the evolution of the local baryon and energy 
densities in the most central cell during the collision of Au-|-Au at P;,eamM= 11-6 GeV/c 
and b=0. Results from the cascade calculations are compared with those from calculations 
using the soft nuclear equation of state with K — 200 MeV and the stiff one with K — 380 
MeV. As expected, the mean field has almost no effect in the early stage of the coUision 
when the kinetic energy is much higher than the potential energy. Significant differences 
appear soon after about 2 fm/c. In the cascade case, it is seen that a maximum baryon 
density of about 9po and a maximum energy density of about 3.6 GeV/fm^ are reached at 
about 4 fm/c. The matter in the high energy density region lasts for about 5 fm/c. The 
volume of the high density region are, however, significantly reduced by the repulsive mean 
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field. Furthermore, the hfetime, maximum baryon and energy densities are reduced to about 
3 fm/c, 7po and 2.6 GeV/fm^, respectively, when using the soft nuclear equation of state. 
With the stiff equation of state the reduction is even larger. Since the stiff equation of state 
violates causality already at about 3po, we will only use the soft equation of state in the 
following. The reduction of the maximum baryon and energy densities due to the mean field 
may be large enough to affect significantly the collision dynamics. 

Additional information about the importance of the nuclear mean field can be obtained 
by comparing the fraction of particles with local energy densities e/ > 2.5 GeV/fm^ from 
calculations with and without the mean field. This is shown in Fig. 22. In the cascade 
calculation about 60% of particles are in the high energy density region at the instant of 
the maximum compression. The soft equation of state reduces it to about 30%, and this 
reduction persists during the expansion phase of the collision. 

C. Thermalization and radial flow 

To what degree the superdense hadronic system is thermalized and what is the time 
scale for achieving thermal equilibrium are among the most important questions in heavy- 
ion collisions at AGS. Answers to these questions are useful both for developing static 
models and for interpreting the experimental data. For example, in thermodynamic models 
one usually assumes that there exists a thermalized freeze-out stage when particles in the 
expanded system cease to strongly interact with each other. In view of the apparent success 
of the thermodynamical approach (e.g. 0]) in explaining a large amount of data from 
heavy-ion collisions at AGS energies we study in this section the question of thermalization. 
In addition, the large thermal and compressional pressure in the superdense hadronic matter 
is found to induce a radial flow during the expansion phase. This is already qualitatively 
clear from our discussions in the previous section, and we study here more quantitatively 
the radial flow. 

We consider all hadrons in a sphere of radius 2 fm around the origin in the center of 
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mass frame for head-on collisions of Au+Au at Pbeam/-^ — 11-6 GeV/c. The degree of 
thermalization can be measured by the ratio 



R = ^'^^"-^/^^^ (99) 



Then i? = 1 is a necessary, although not a sufficient, condition for thermal equilibrium. It 
is also a measure of the stopping power. In the bottom window of Fig. 23, we show this 
ratio as a function of time. For all three calculations based on the cascade, the soft and the 
stiff equation of state, the ratio reaches one at about 10 fm/c. From our early discussions 
on the reaction rates for various channels, it is clear that this is about the time when most 
hadronic coUisions have stopped. It is therefore reasonable to use the concept of thermal 
equilibrium at the time of freeze-out. 

It is interesting to note that at about 4 fm/c when the system reaches its maximum 
compression the ratio R is about 0.7, indicating that the superdense hadronic system has 
not reached thermal equilibrium. To estimate the temperature of particles in the above 
sphere at the time of thermalization, we show in the middle window the scaled average 
kinetic energy per particle < >. At freeze-out, this quantity can give a good estimate 
of the temperature. Before freeze-out, the evolution of this quantity reflects the coUision 
dynamics. It is seen that the kinetic energy decreases quickly in the compression phase due 
to the creation of baryon resonances and other particles as well as various scattering. In the 
expansion phase, high energy particles near the surface of the sphere move out continuously 
due to the radial flow. The average kinetic energy in the sphere thus gradually decreases in 
the expanse phase. At the time of freeze-out this quantity is about 120 MeV and shghtly 
depends on the model. This value is very close to the temperature extracted from the inverse 
slope of the measured particle spectra. It is interesting to note that both the time scale for 
thermalization and the temperature at freeze-out are not very sensitive to the equation of 
state. However, the temperature from the cascade is about 30 MeV higher than that from 
calculations with the mean fleld due to the neglect of the potential energy in the former. 

In our calculations, the system has undergone sufficient expansion before reaching ther- 
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mal equilibrium at freeze-out, and this is in agreement with the conclusion reached from 



studies of the AGS data based on the thermodynamical model ||5^. More quantitative in- 
formation about the expansion can be seen from the upper window where the number of 
particles in the central sphere is plotted as a function of time. General features of this 
plot are the same as that of the central density plot discussed earlier. It is seen that the 
number of particles in the sphere at freeze-out is only about 30% of that at the maximum 
compression indicating a strong expansion before the freeze-out. 

The velocity profile of the radial flow at freeze-out can be expressed by 

0(r) = Jfj:'^ ■ f, (100) 

where A^,- is the number of particles in a spherical shell of size Ar = 0.5 fm at a radius r. 
The momentum pi of each particle has a flow component and a random thermal component 

Pi = Pif + Pit- (101) 

In the limit that there are a large number of particles in each shell, the summation in Eq. 
( |100| ) cancels out the thermal part, i.e. 

Nr=O0 -> 

^.11 = 0, (102) 

i=l ''^i 

and Eq. ( [I00| ) reduces thus to the flow velocity Pf{r), 

/3(r)^/3/(r) = ''lf (103) 
i=i ^« ''"^ 

Since we usually use a large number of test particles, Eq. (|100|) gives a good estimate of 
the radial flow profile. In Fig. 24, we show the estimated flow proflle and the corresponding 
radial density distributions for baryons and pions. The left window is for the cascade and 
the right one is for calculations using the soft equation of state. The radial flow velocity 
is proportional to the radius, but there seem to be two distinct slopes at radii smaller and 
larger than about 6 fm. Another interesting feature is that pions seem to flow a little faster 
than baryons due to their lighter masses. It is also seen that the difference between the 

34 



flow profiles from calculations with and without the mean field is rather small, which may 
indicate that the thermal pressure dominates at freeze-out. As seen from the upper window, 
more particles are at larger radii due to the repulsive mean field, so more energy are stored 
in the radial flow in the case with the mean field. 

For comparisons, we mention here that results of a recent hydrodynamical model analysis 
of a large amount of data from central collisions of Si+Au at Pbeam/^ = 14.6 GeV/c are 
consistent with the above discussions. Assuming a freeze-out radius of about 7 fm, a freeze- 
out temperature of 120 (140) MeV and a linear flow velocity profile, the extracted maximum 



flow velocity is about 0.58 (0.5) c [54|, which agrees with the flow profile we have extracted. 



D. Stopping power and inclusive observables 

The formation of a quark-gluon plasma or a mixed phase of hadrons, quarks and gluons 
is expected to reduce the pressure of the system and leads to a softened equation of state. 
It is therefore interesting to search for experimental observables that are sensitive to the 
nuclear equation of state. In this section, we study the sensitivity of inclusive, single particle 
observables to the equation of state by comparing calculations with and without the mean 
field. However, we should mention again that present models for heavy-ion collisions at AGS 
energies are probably inconsistent. Namely, we have shown in the above that the conditions 
for forming the quark-gluon plasma may have been fulfilled in central collisions at AGS 
energies. Without actually including the phase transition to the quark-gluon plasma and 
the hadronization back to the hadronic phase in the model, we are limited to discuss only 
the hadronic physics. 

Specifically, we shall study the rapidity and transverse momentum distributions of pro- 
tons and pions. We leave the study of kaon distributions to the next section in central 
collisions of Au+Au at pbeam/A =11.6 GeV/c. For references, we also compare with prelim- 
inary data from the E802/E866 collaboration |^. These data were taken using the kinetic 
energy in the zero degree calorimeter as a trigger on the impact parameter. The data for 
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central collisions roughly correspond to the upper 4% of the total cross section or impact 
parameters smaller than about 3 fm. Since the data still have a 15% systematic errors be- 
sides the statistical ones, we have not let the generated events going through the zero degree 
calorimeter to simulate the experimental trigger conditions. Instead, we perform calcula- 
tions with impact parameters smaller than 3 fm and compare the data with calculations 
using the impact parameter of 2 fm which is about the average impact parameter of the 
central collisions. This is probably sufficient for our current purposes. 

First, let us study the proton transverse momentum distribution. In Fig. 25, we show 
the calculated spectrum by the round symbols in the left window using the cascade and 
the right window using the soft equation of state. The least square fit to the data of ref. 



55| are shown by the solid lines. In the cascade case, the calculation agrees with the data 
rather well except in the low transverse momentum region. The larger theoretical values 
are probably due to the fact that we count all protons in the final state although some 
of them are still bound in light clusters, while the data are for protons only. To estimate 
the number of bound protons, we apply a cut to the local baryon density such that only 
protons with pi < 0.2po are counted. The spectra calculated in this way are shown with the 
cross symbols which agree well with the data even at low transverse momenta. The spectra 
calculated with the soft equation of state show an increase of about 25% mainly at large 
transverse momenta compared with that from the cascade as a result of the repulsive mean 
field at high densities. This difference, however, is comparable to the combined statistical 
and systematical errors of the data. Nevertheless, it is promising that the proton spectra 
at high transverse momenta may be useful for studying the equation of state at the level of 
sensitivity of about 25%. As we will show later, a much clearer sensitivity at the level of 
250% can be obtained from an event by event analysis of the baryon transverse collective 
fiow. 

In Fig. 26, we show the calculated proton rapidity distributions. As one would expect, the 
mean field leads to less stopping. This observation is consistent with our early discussions 
about the mean field effects on the creation of superdense hadronic matter. It is seen 
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that without applying the local density cut at 0.2po, the mean field reduces the height but 
increases the width of dN/dy at midrapidity by about 15%. With the local density cut, 
the calculation with the mean field can better reproduce the data, but again the difference 
between the two calculations is comparable to current systematic errors in the data. As 
we have stated earlier, the above calculations are carried out at a fixed impact parameter 
of 2 fm. To quantify the impact parameter dependence of the stopping power in central 
collisions, we show in Fig. 27 the proton rapidity distributions at impact parameters of 1 fm 
and 3 fm. In both calculations with and without the mean field, the height of the rapidity 
distribution at the midrapidity increases by about 25% as the impact parameter changes 
from 3 fm to 1 fm. With an impact parameter of 3 fm, a clear bump due to spectators at 
the projectile and target rapidities can be seen. The variation of the rapidity distribution 
at these two impact parameters caused by the mean field is again at the level of 15%. 

Pions are copiously produced in heavy-ion collisions at AGS energies and are expected to 
play an important role in the stopping and thermal equilibration of the system. The study 
of pion rapidity and transverse momentum spectra provides thus useful information about 
the collision dynamics, the properties of the superdense hadronic matter, and the particle 
production mechanisms. In Fig. 28 and Fig. 29, we compare calculated tt"^ and 7r~ transverse 
momentum spectra using the cascade and soft equation of state with each other and also with 
the least square fit to the data (solid lines) . The rapidity range for each spectrum is indicated 
in the figure. It is seen that both calculations can reproduce the data reasonably well except 
at rapidities far from the center of mass rapidity, which is about 1.72 for the Au+Au reaction 
at Pbeam/^=ll-6 GcV/c. Again, the variation between calculations using the cascade and the 
soft equation of state is on the order of 20% which is comparable to current experimental 
errors. The statistics in the calculations can be significantly improved by mixing events 
from different impact parameters. But main features of the spectra are expected to remain 
unchanged. We will leave a more detailed study on pion spectra, including the contribution 
of resonance decays at freeze-out and the Coulomb effects on the difference between tt"*" and 
7r~ spectra to a separate paper. 
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In Fig. 30, we compare the rapidity distributions of 7r+ calculated at an impact parameter 
of 2 fm using the cascade and the soft equation of state with the experimental data. The 
calculation with the soft equation of state seems to be closer to the data, but there exists 
some discrepancies. It is worth mentioning that the experimental rapidity distribution was 
obtained by first fitting the transverse momentum spectra with an exponential function or 
a sum of two exponential functions and then extrapolating and integrating over the whole 
rapidity range. While in our calculations we simply count pions in each rapidity bin without 
any possible bias from the extrapolation process. It is seen that the difference between 
calculations with and without the mean field is again on the order of 20%, which is rather 



similar to the situation at the Bevalac or SIS/GSI energies |1T7| , |43[| . The impact parameter 
dependence of the pion rapidity distribution is shown in Fig. 31. The dependence here is 
rather weak compared to that of the proton rapidity distribution. This is mainly due to 
the fact that there is an almost complete overlapping between the target and the project in 
collisions with impact parameters smaller than 3 fm and the large number of pions produced 
in these collisions. 



E. Mechanisms for production 

Mechanisms for kaon production at AGS energies have been a subject of intensive studies 
in the past few years. In particular, the enhancement of /tt^ ratio in heavy-ion collisions 
with respect to pp and pA reactions has been attracting much attention due to the possible 
role of strangeness as a signature of the quark-gluon plasma fS^JS^. However, no definite 
conclusion has been reached due to the complexity of the collision dynamics at AGS energies 
and our incomplete knowledge about the properties of hadrons and their interactions in hot 
dense matter. 

In Fig. 32, we show the contributions to production from various collision channels 
as a function of the center-of-mass energy of the colliding pairs in Au+Au collisions at 
Pbeam/A = 11.6 GcV/c. It is sccu that large contributions come from meson+baryon and 
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meson+meson collisions, besides that from baryon+baryon collisions. This is in agreement 



with the predictions from ref. |^ based on the hydrochemical model calculations. More 
quantitatively, about 35%, 40% and 25% of the produced kaons come from baryon+baryon, 
meson+baryon and meson+meson collisions, respectively. The relative importance of dif- 
ferent K'^ sources seen here does not agree completely with results from either the RQMD 
or the ARC. Since there are many assumptions in each model, it is not apparent to us at 
this time what is the cause for the difference among these models. It is certainly of great 
interest to resolve this difference in future studies. However, this can only be possible when 
all assumptions used in all models are made clear. 

For ease of future comparisons with other models and also for studying the kaon pro- 
duction dynamics, we show in Fig. 33 the production rate as a function of time for 
various channels. The calculation is carried out for head-on collisions at Pbeam/^ = 11-6 
GeV/c using the cascade. It is seen that at t ~ 1.0 fm/c the production rates due to nu- 
cleon+nucleon, meson+baryon (including resonances) and nucleon+baryon resonance are 
almost the same. Soon after that the production rate is dominated by meson+baryon 
collisions which also last a longer time. The production rates due to nucleon+baryon 
resonance, nucleon+nucleon and baryon resonance+baryon resonance decrease to zero after 
about 6 fm/c. Meson+meson collisions come in later, but they contribute significantly to 
the production especially during the expansion phase of the collision. The time depen- 
dence of production studied here may be useful for further studies on the production 
mechanisms for kaons. Moreover, kaons produced at different times through different chan- 
nels may also contribute to different parts of the kaon spectra. We, however, leave the more 
detailed discussions on this subject to future studies. 

We now go back to the question of comparing kaon spectra calculated with and without 
the mean field with the preliminary data from the E802/E866 collaboration. Fig. 34 shows 
the comparison on the kaon transverse momentum distributions. The solid lines in both win- 
dows are the exponential fit to the data to guide the eye. It is seen that both calculations 
agree with the data reasonably well. However, there are some discrepancies between cal- 
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culations and experimental data especially at high transverse momenta. Within error bars 
of the calculations the two calculated spectra are almost indistinguishable. Fig. 35 shows 
the rapidity distributions. The data at rapidities larger than 1.72 are the reflections 
of those at rapidities less than 1.72. Except at rapidities near the center-of-mass rapidity 
the two calculations are very close to each other. Even at the midrapidity the sensitivity 
to the equation of state is only about 15%. This situation is rather different from that at 
SIS/GSI energies where subthreshold kaon production is very sensitive to the nuclear equa- 



tion of state ||58|-p4l. Since kaon production at AGS energies is high above the threshold, 
the sensitivity of kaons to the equation of state is not expected to be very different from that 
of nucleons and pions. To be complete, we show the impact parameter dependence of kaon 
production in Fig. 36. The sensitivity of kaon multiplicity to the impact parameter is at the 
level of 20% in central collisions, which is rather similar to that of pions. The same level 
of sensitivity to both the impact parameter and the equation of state in central collisions 
makes it even harder to learn about the equation of state from these observables. In contrast 
to the inclusive observables discussed above, collective observables are very sensitive to the 
equation of state as we shall show in the next section. 

F. Baryon transverse flow 

Hinted by the flndings at Bevelac energies that the collective variables or correlation 
functions, unlike single particle observables, are very sensitive to the nuclear equation of 



state, we now turn to the analysis of the baryon transverse collective flow |^3[. Here we use 



the standard flow analysis in transport models as in ref. [|TT|. It has been shown recently in 



ref. [0 that an improved flow analysis with the explicit conservation of reaction plane and 



angular momentum in an individual hadron-hadron collision results in an increase of the 
transverse momentum by about 8% to 23% at Bevalac energies. This small enhancement is 
not expected to affect our discussions and conclusions. 

First, we show in Fig. 37 the average transverse momentum of nucleons in the reaction 
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plane as a function of rapidity for the collision of Au+Au at Pbeam/^ =11.6 GeV/c and at 
impact parameters of 2, 6 and 10 fm. It is seen that the collective flow is the strongest 
in midcentral collisions in both calculations. Most importantly, significant differences exist 
between calculations with and without the mean field for collisions at all three impact 
parameters. In particular, the flow parameter defined as 



is about a factor of 2.5 larger in the case with the mean field. The strength of the so-called 
"bounce-off" effect at target or projectile rapidities is also much stronger in calculations with 
the mean field. To measure the strength of the "bounce-off" effect, we define the average 
total in-plane transverse momentum as 



To see how the collective flow is a sensitive probe of the reaction dynamics in the high 
density region, we show in Fig. 38 and Fig. 39 the flow parameter F and the total in- 
plane transverse momentum < P-^ > as functions of time for Au-|-Au reactions at impact 
parameters of 2 fm and 6 fm, respectively. For comparisons, the central densities are also 
shown in the lower windows of Fig. 38 and Fig. 39. It is clearly seen that the flow is mainly 
generated in the high density region and does not change in the expansion phase. The 
increase of the total in-plane transverse momentum during the expansion phase is much 
more substantial in the case with the mean fleld than in the cascade case. This indicates 
that the mean field effect is large in the expansion phase as we have expected. The total 
in-plane transverse momentum reaches its final value after about 15 fm/c. The ratio of final 
flow parameters in the calculations with and without the mean field is about 2.5 at both 
impact parameters. This difference is large enough for distinguishing the predictions from 
different models. For collision at an impact parameter of 6 fm, the fiow parameter decreases 
slightly before reaching its final value. This is mainly due to the refiection of hot baryons 
from the cold spectator nucleons in collisions at large impact parameters. 




(104) 




(105) 
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It is interesting to note that a clear signature of transverse collective flow has been 
discovered recently in Au+Au reactions at Pbeam/ ^= 11-4 GeV/c by the E877 collaboration 
p5| . In this experiment the distribution of the normalized transverse energy dipole moments 
is found to shift systematically toward flnite values for more central collisions. We are 
planning to make a detailed comparison with these data in the near future. Furthermore, 
the preliminary data of the E877 collaboration indicates a smaller transverse momentum in 
the reaction plane than the predictions from both present and other models ||66| , |67[| , which 
may indeed indicate that the equation of state is significantly softened due to the formation 
of the quark-gluon plasma. 



IV. SUMMARY 

In summary, we have developed a new relativistic transport model (ART 1.0) for heavy- 
ion collisions at the AGS energies. This model is an extension of the very successful BUU 
model for intermediate energy heavy-ion collisions by including more baryon and meson res- 
onances as well as interactions among them. Comparing to the well-known ARC model [Q, 
the present model has more resonances which are found to be important for the collision dy- 
namics and particle production at AGS energies. Moreover, we have used consistently finite 
mass distributions for these resonances which have significant effects on particle spectra. In 
addition, a selfconsistent mean field for baryons is included. However, comparing to another 
well-known RQMD model 0, where all resonances up to mass of 2.0 GeV are explicitly 
included, we do not have as many baryon resonances. Nevertheless, we have partially in- 
cluded the effects of these resonances by using meson+baryon cross sections calculated from 
the implicit formation of higher resonances. As we are mainly interested in developing a 
reliable model for energies up to about 15 GeV/nucleon, we do not invoke the string/rope 
mechanisms for particle production as in RQMD. 

Due to the lack of some hadron-hadron collision data, the input to the model still has 
some uncertainties. Several assumptions used in the model may need to be improved in 
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future versions of the model. In particular, interactions among resonances are largely un- 
known and have been assumed to be the same as nucleon-nucleon interactions. Additional 
hadrons, such as p, A, etc, are to be included in the model. At AGS energies the resonance 
population is as large as that of nucleons in the compression phase, the assumption of the 
same mean field for nucleons and resonances may need to be improved. Currently, the finite 
formation time for produced particles has only been partially taken into account through 
the finite lifetime of the resonances. A more consistent method may be necessary. Also, it is 



of interest to include the mean field potential for pions |l68| , |69[| and kaons |]70[ as well as the 
Bose enhancement factor for pions [^. The most inconsistent, but also probably the most 
interesting feature shared by this model and others is the prediction that the conditions for 
forming the quark-gluon plasma have already been reached in heavy-ion collisions at AGS 
energies. In particular, at the maximum baryon density of about 9po and meson density of 
4.5po hadrons are squeezed on top of each other and probably melt into the quark-gluon 
plasma. At this high density hadronic models naturally breaks down. It remains a great 
challenge to actually include the quark and gluon degrees of freedom into the model and to 
study the phase transition. 

With the above reservations about the validity of the model, we have studied several 
aspects and issues of heavy-ion reactions at AGS energies. In particular, we have explored 
effects of the nuclear mean field on the collision dynamics, the formation of superdense 
hadronic matter, single particle observables as well as the radial and transverse collective 
flow. It is found that the mean field significantly reduces the maximum energy and baryon 
densities of the superdense hadronic matter formed in the collision. Inclusive single particle 
observables are sensitive to the equation of state only at the level of 20%. On the other hand, 
the transverse collective flow of baryons are found to be very sensitive to the equation of 
state. In particular, including a soft nuclear equation of state increases the flow parameter by 
a factor of 2.5 which is large enough to distinguish the predictions from different theoretical 
models. Moreover, this sensitivity also indicates that the collective flow can be a useful probe 
of the formation of the quark-gluon plasma at AGS energies. Preliminary data have shown 
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a smaller baryon transverse flow in the reaction plane than model predictions, and this may 
indeed indicate the significant softening of the equation of state due to the formation of the 
quark-gluon plasma. It is very encouraging that the present model has been very useful for 
studying the interesting physics in high energy heavy-ion collisions. 
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TABLES 



TABLE I. Isospin cross section parameters 



parameter an aio ctqi 

a 3.772 15.28 146.3 

/3 1.262 

rno(McV) 1188 1245 1472 

r(MeV) 99.02 137.4 26.49 



TABLE IL Higher A^* resonances 



R 


JP 


Mass 


Width 


Ntt 


Att 


Np 


^13 


3- 
2 


1.520 


0.125 


0.550 


0.250 


0.200 


Su 


1- 
2 


1.650 


0.150 


0.600 


0.050 


0.180 




5- 
2 


1.675 


0.155 


0.380 


0.580 


0.030 




5 + 
2 


1.680 


0.125 


0.600 


0.125 


0.125 




3- 
2 


1.700 


0.100 


0.100 


0.380 


0.100 


Pu 


1 + 
2 


1.710 


0.110 


0.150 


0.100 


0.200 


Pl3 


3 + 
2 


1.720 


0.200 


0.150 


0.100 


0.530 


Fn 


7 + 
2 


1.990 


0.290 


0.050 


0.060 


0.340 
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TABLE III. Higher A resonances 



R 


JP 


Mass 


TTT* 1j_1 

Width 


A r 

A'TT 




Np 


-F33 


3 + 
2 


1.600 


0.250 


0.350 


0.450 


0.050 


S31 


1 — 
2 


1.620 


0.160 


0.300 


0.600 


0.070 


D33 


3 ~ 
2 


1.700 


0.280 


0.150 


0.698 


0.150 


S31 


1 — 
2 


1.900 


0.150 


0.100 


0.050 


0.450 


P35 


5 + 
2 


1.905 


0.300 


0.100 


0.250 


0.450 


^31 


1 + 
2 


i.yiu 


U.zzU 


n oon 
U.zZU 


U.Uoy 


U.UoU 


P33 


3 + 
2 


1.920 


0.250 


0.200 


0.190 


0.080 


D35 


5- 
2 


1.930 


0.250 


0.090 


0.200 


0.120 


F37 


7 + 
2 


1.950 


0.240 


0.400 


0.130 


0.080 
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FIGURE CAPTIONS 



Fig. 1 pp inelastic cross sections as functions of center-of-mass energy ^/s. 

Fig. 2 pp — > ppcu and pp — > ppp° cross sections as functions of center-of-mass energy. 

Fig. 3 The + p elastic cross section as a function of center-of-mass energy. The solid 
line is the experimental data and the dotted line is the contribution from the formation 

of A, iV*(1440) and iV*(1535) resonances. The dashed line is the difference between 
the solid and dotted lines and is attributed to direct tt -|- scattering. 

Fig. 4 7r~^ + p elastic cross section as a function of center-of-mass energy. The solid line 
is the experimental data and the dotted line is the contribution from the formation of 
the A resonance. The dashed hne is the difference between the solid and dotted lines 
and is attributed to direct tt + N scattering. 

Fig. 5 The +p{n) elastic cross section as a function of center-of-mass energy calculated 
from the resonance model. 

Fig. 6 The tt"*" -|- A~(A°) elastic cross section as a function of center-of-mass energy 
calculated from the resonance model. 

Fig. 7 The n'^ + p inelastic cross section as a function of the center-of-mass energy ^/s. 

Fig. 8 The isospin-averaged tttt elastic (upper window) and inelastic (lower window) cross 
sections as functions of center-of-mass energy y/s. 

Fig. 9 The angular distribution of NN — > A'"A reaction at different center-of-mass 
energies. 

Fig. 10 Transverse momentum distributions for protons (upper window) and pions (lower 
window) from pp collisions at P6eam=15 GeV/c and 3 GeV/c. Solid lines are the scaling 
results valid at high energies as discussed in the text. 
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Fig. 11 The average energy per nucleon (upper window) and the adiabatic sound velocity 
(lower window) as functions of density in nuclear matter at zero temperature. 

Fig. 12 Reaction rates for various channels in nucleon+nucleon coUisions as functions of 
time in head-on collisions of Au+Au at Pbeam/^ = 11-6 GeV/c. 

Fig. 13 Reaction rates of nucleon+baryon resonance, baryon resonance+baryon resonance 
and meson+meson collisions as functions of time in head-on collisions of Au-|-Au at 

Pbeam/A = 11.6 GcV/c. 

Fig. 14 Rates of baryon resonance decays and formations in meson-|-baryon coUisions as 
functions of time in head-on collisions of Au-|-Au at Pbeam/^ — 11-6 GeV/c. 

Fig. 15 The evolution of the tt, A, and N* multiphcities in head-on coUisions of Au-|-Au 
at Pbeam/A =11.6 GeV/c. 

Fig. 16 Collision number distributions as functions of center-of-mass energy in head-on 
collisions of Au-|-Au at Pbeam/^ =11.6 GeV/c. 

Fig. 17 The evolution of local baryon density in the reaction plane in head-on collisions 
of Au-fAu at Pbeam/A = 11.6 GcV/c. 

Fig. 18 Same as Fig. 17 for the scalar baryon density. 

Fig. 19 Same as Fig. 17 for the local meson density. 

Fig. 20 Same as Fig. 17 for the local energy density. 

Fig. 21 The evolution of central local baryon and energy densities in the collision of 
Au-|-Au at Pbeam/^ =11.6 GcV/c and b=0 calculated with the cascade, the soft, and 
the stiff equation of state. 

Fig. 22 The evolution of the fraction of particles with local energy density higher than 
2.5 GeV/fm^ in head-on collisions of Au-|-Au at P^eamM =11-6 GeV/c. 
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Fig. 23 (Upper window): the evolution of the particle number in a sphere with a radius 
of 2 fm around the center of mass. (Middle window): the evolution of the average 
kinetic energy < ^Ek > per particle in the sphere. (Lower window): the evolution of 
the stopping ratio R in head-on collisions of Au+Au at PbeomM =11-6 GeV/c. 

Fig. 24 (Upper windows): Radial density distributions for baryons and pions at freeze- 
out calculated using the cascade (left) and the soft equation of state (right). (Lower 
windows): the radial flow velocity profiles at freeze-out in head-on collisions of Au-(-Au 

at Pbeam/A ^11.6 GcV/c. 

Fig. 25 Transverse momentum distributions of protons in the collision of Au-|-Au at 
Pbeam/A =11.6 GcV/c and at an impact parameter of 2 fm. The solid line is the least 
square fit to prehminary data from the E802 collaboration. See discussions in the text. 

Fig. 26 Rapidity distributions of protons in the coUision of Au-|-Au at Pbeam/A =11.6 
GeV/c and at an impact parameter of 2 fm. 

Fig. 27 Comparison of rapidity distributions of protons in the collision of Au-|-Au at 
Pbeam/A =11.6 GcV/c and at impact parameters of 1 and 3 fm. 

Fig. 28 Transverse momentum distributions of pions from the cascade model for Au-|-Au 
collisions at Pbeam/A =11.6 GeV/c and at an impact parameter of 2 fm. Solid lines 
are least square fits to preliminary data from the E802/E866 collaboration. 

Fig. 29 Same as Fig. 28 for transverse momentum distributions of pions. 

Fig. 30 Rapidity distributions of 7r+ in the collision of Au-|-Au at Pbeam/A =11.6 GeV/c 
at an impact parameter of 2 fm. 

Fig. 31 Comparison of rapidity distributions of 7r+ in the coUision of Au-|-Au at 
Pbeam/A =11.6 GcV/c and at impact parameters of 1 and 3 fm. 
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Fig. 32 Energy distributions of K'^ producing collisions as functions of center-of-mass 
energy in head-on collisions of Au+Au at Pbeam/^ =11.6 GeV/c. 

Fig. 33 The evolution of K'^ production rates in head-on collisions of Au-|-Au at Pbeam/^ — 
11.6 GeV/c. 

Fig. 34 Transverse momentum distributions of K'^ in the collision of Au-|-Au at 
Pbeam/A =11.6 GcV/c at an impact parameter of 2 fm. Solid lines are least square fits 
to prehminary data (solid symbols) from the E802/E866 collaboration. 

Fig. 35 Rapidity distributions of K'^ in the coUision of Au-|-Au at Pbeam/A =11.6 GeV/c 
at an impact parameter of 2 fm. 

Fig. 36 Comparison of rapidity distributions of K'^ in the coUision of Au-|-Au at 
Pbeam/A =11.6 GeV/c and at impact parameters of 1 and 3 fm. 

Fig. 37 Baryon average transverse velocity in the reaction plane as a function of rapidity 
for Au-|-Au collisions at P6eam/^=ll-6 GeV/c and at impact parameters of 2, 6, and 
10 fm. 

Fig. 38 Time evolution of the flow parameter, the central density and the total in-plane 
transverse momentum in the collision of Au-|-Au at P6eam/^=ll-6 GeV/c and at an 
impact parameter of 2 fm. 

Fig. 39 Same as Fig. 38 at an impact parameter of 6 fm. 
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